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Abstract. Theoretical progress on the for¬ 
mulation and classification of the quark 
and gluon orbital angular momenta (OAM) 
is reviewed. Their relation to parton dis¬ 
tributions and open questions and puz¬ 
zles are discussed. We give a status report 
on the lattice calculation of the parton ki¬ 
netic and canonical OAM and point out 
several strategies to calculate the quark 
and gluon canonical OAM on the lattice. 

PACS. 12.38.-t Quantum chromodynam¬ 
ics - 14.20.Dh Protons and neutrons - 
12.38.Gc Lattice QCD calculations 

1 Introduction 

One of the key questions in hadron physics is to under¬ 
stand how the spin of the proton originates from the par- 
ton spin and orbital motion. Based on the naive quark 
model picture, supported among others by its success in 
hadron spectroscopy, it was expected that the quark spin 
contribution should account for mosiQof the proton spin. 
Deep-inelastic scattering experiments have been conducted 
and obtained the surprising result that the quark spin 
contribution is actually quite small, about 25% [Tj. Rela¬ 
tivistic quark models argue that the quark orbital motion 
reduces significantly the quark spin contribution, but in 
practice this effect falls short in explaining such a small 
value. According to Quantum ChromoDynamics (QCD), 
gluon polarization also contributes to the proton spin bud¬ 
get. This contribution appears to be of the order of 35% [351 

1 In relativistic quark models, the quark orbital motion can 
contribute significantly to the proton spin budget, but the 
quark spin remains by far the dominant contribution m- 


I139j . leaving about 40% which should come from the or¬ 
bital angular momentum (OAM) of quarks and gluons. 
This number is large and reflects the relativistic nature of 
the quark-gluon bound state. It also stresses the paramount 
importance of actually measuring the OAM contribution. 

The parton OAM is more difficult to access than the 
parton helicity owing to the fact that it measures a corre¬ 
lation between position and momentum, and therefore re¬ 
quires some information about parton phase space. More¬ 
over, it also depends on how the total AM is divided into 
separate quark and gluon contributions which is intrinsi¬ 
cally ambiguous due to quark-gluon couplings. It has long 
been thought that only the so-called kinetic decomposi¬ 
tion of the proton spin makes sense because its quark and 
gluon contributions can be extracted from experiments 
and computed on the lattice. Recent theoretical and ex¬ 
perimental progress have however shown that the canoni¬ 
cal decomposition can also be accessed experimentally and 
computed on the lattice, though in a more complicated 
way. Kinetic and canonical decomposition are two com¬ 
plementary decompositions of the nucleon spin with their 
own advantages and disadvantages. For more detailed dis¬ 
cussions about the angular momentum (AM) decomposi¬ 
tion in gauge theories, see the recent reviews [94lll56j . For 
a theoretical and experimental overview of optical AM, 
see [20] . 

The aim of the present paper is to summarize criti¬ 
cally the current theoretical knowledge about the parton 
OAM and give some prospects for future interesting de¬ 
velopments. The paper is organized as follow. We quickly 
review in sec. [2] the recent theoretical development about 
the proton spin decomposition and discuss both kinetic 
and canonical decompositions. In Sec. [3j we show how 
the different contributions, and in particular the parton 
OAM, can be expressed in terms of parton distributions 
and therefore how they can be accessed experimentally. 
We then discuss in detail in Sec. [4] several issues and puz¬ 
zles related to the definition and calculation of OAM. We 
also comment critically various claims found in the litera¬ 
ture. In Sec.[5j we present in a nutshell some lattice calcu¬ 
lations of the various AM contributions, starting first with 
the Mellin moment approach followed then by the more 
recent large momentum approach. We suggest in Sec. [6] 
that a realistic effective theory of baryons consistent with 
chiral symmetry could provide valuable quantitative in¬ 
sights concerning the question of the proton spin decom¬ 
position. In Sec. [7] we present the concept of spin-orbit 
correlation and show that it is closely related to parton 
OAM. In particular, all the discussions and developments 
regarding OAM find a natural counterpart when studying 
the spin-orbit correlation. Finally, we summarize in Sec. [8] 

2 The gauge-invariant angular momentum 
decomposition 

Generically, the total AM operator in QCD can be decom¬ 
posed into five gauge-invariant contributions [T54j 

J = S^ + S G + Ll in + Lg n + L pot , 


(1) 
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where 

S q = J A z r^\Exl), 

S G = J d 3 4 rE a x Ap hys , 

Kin = J d 3 ripHr x iD)ip, (2) 

Lg a = ~ j d 3 rE a \r x V a p l e )A% ys , 

Lpot = - J d 3 rp a r x A“ hys 

are the quark spin, gluon spin, quark kinetic]^] OAM, gluon 
canonical OAM and potential OAM, respectively. 

Following Chen et al. nm, the gauge potential has 
been decomposed into a pure-gauge Held and a “physical” 
field 

^(r) = A£ ure (r) + A£ hys (r) (3) 

such that 

A v — c) v —in\A^ A v 1=0 (A) 

J pure w -^pure w ' r:L pure Impure’ -'■‘pureJ V^v 

Under a gauge transformation, A^ ure is assumed to trans¬ 
form as a connection, i.e. like the gauge potential A 

A pnre( r ) ^ d( r ) Impure O’) + (5) 

It follows that Ap hys , which is the difference between two 
connections, transforms homogeneously 

^physM h + 9 (r)A» hys (r)g-\r). (6) 


the potential OAM is combined with the quark kinetic 
OAM, one obtains the quark canonical OAM 


r q — r q 

^can -^kin 


-'pot 


/dW(r- x iD pure )ip, 


(9) 


leading to the so-called gauge-invariant canonical or Chen 
et al. decomposition fTTo~llTT?T| 

J = S q + S G + L q can + Lg n (10) 

which can be seen as a gauge-invariant version 3 ] (or ex¬ 
tension) of the Jaffe-Manohar decomposition [72]. When 
the potential OAM is combined with the gluon canonical 
OAM, one obtains the gluon kinetic OAM 

r G _ tG i r 

-^kin — -‘-'can T -‘-'pot 

= J d 3 rrx[(A; hys xS;S re )x£ k ] 

with D£ ure = \(V^+V£ me ) m, leading to the so-called 
gauge-invariant kinetic or Wakamatsu decomposition (1541 

J = S q + S G + Ll in + Lg n (12) 

which generalizes the Ji decomposition m 

J = S q + Ll in + Jg n , (13) 

where 

Jfin = + L G n = j d 3 r [r x (E a x B a )\ (14) 


involves E a x B a , the Poynting vector of the glue. 


The pure-gauge covariant derivatives are defined like the 
ordinary covariant derivatives with the gauge potential A^ 
replaced by the pure-gauge field A£ ure 

-Dpure = ^ - igA pi]re , 

V£ uie = d»-ig[A£ uie , }. U 

This ensures the explicit gauge invariance of the expres¬ 
sions in Eq. ([2]). It turns out that the Chen et al. ap¬ 
proach is just a particular case of the standard background 
field method with the background field A puro chosen to be 
flat 111 2 111131 . i.e. satisfying the pure-gauge condition in 
Eq. (jijb 

Because of the QCD equations of motion 

p a = g^t a ip = T> ab ■ E b , (8) 

the interpretation of the potential OAM as a quark or 
gluon contribution is ambiguous. It represents in some 
sense the change in OAM as the quark moves through 
the color field created by the spectator system |32j. When 

2 The kinetic OAM is sometimes called mechanical or dy¬ 
namical OAM in the literature. 


The gauge invariance of the terms in Eq. Q seems to 
contradict the textbook claim that a gauge-invariant de¬ 
composition of the gluon total AM into spin and OAM 
contributions is impossible [171173] . The contradiction is 
only apparent though since textbooks deal actually with 
local fields, whereas A(j ure and Ap hys are intrinsically non- 
locag finimT] . As a result, manifest Lorentz covariance 
is lost. Moreover, the pure-gauge condition Q is not suf¬ 
ficient for fixing uniquely the form of A(( ure and Ap hys , 
leading to a new freedom at the level of the QCD La- 
grangian 

-dpureO) ^ A purei r ) + ^C(r), 

< hys (r)^ A ; hys (r)-d^C(r), 

referred to as the Stueckelberg transformation [145111101 
[45] and related to the concept of background dependence 

3 By gauge-invariant extension, we mean a gauge-invariant 
quantity which gives the same physical results as a gauge- 
variant one considered in a certain fixed gauge. 

4 The fields A£ ure and A£ hya can be expressed solely in terms 
of the full gauge potential A M , but this requires the knowledge 
of A^ at least over a finite neighborhood. 
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|113j . In theory, this freedom is not welcome as it makes 
the initial decomposition intrinsically ambiguous. In 
practice, however, the experimental conditions usually pro¬ 
vide us with a natural choice for the background field 
pi ire mmm^ giving e.g. a preferred spatial direc¬ 
tion and explaining why the photon spin and OAM are 
routinely measured in opticej^] [1531IT91I20] . In high-energy 
experiments where the proton internal structure is probed, 
the off-shell probe provides the natural direction [!37j. In 
other words, even if the fundamental Lagrangian is Stueck- 
elberg invariant (i.e. background independent), physical 
quantities extracted within a controlled expansion scheme 
are often associated with a natural fixed background. 

In summary, a complete gauge-invariant decomposi¬ 
tion of the proton spin makes sense and is feasible, but 
its precise form depends on how one probes the system. 
There is also an intrinsic ambiguity in decomposing the 
total OAM into quark and gluon contributions, leading 
to essentially two types of decompositions: canonical and 
kinetic. This ambiguity can however be avoided by refrain¬ 
ing from interpreting the potential OAM as either quark 
or gluon contribution. 


3 Relations to parton distributions 


Using the Belinfante-Rosenfeld version of the QCD energy- 
momentum tensor, Ji m derived a remarkable relation 
between the quark/gluon total kinetic AA0 and twist-2 
generalized parton distributions (GPDs) 

n) = §1 dx x[H q,G (x, 0,0) + E q,G (x, 0,0)]. (16) 


This relation has been obtained for the longitudinal com¬ 
ponent of the total AM Jl = J ■ P/\P\ and does not de¬ 
pend on the magnitude of the proton momentum ipi m 
1941 . The corresponding relation for the transverse compo¬ 
nent of the total AM Jt has been obtained by Leader [92] 
and is momentum dependent, a consequence of the fact 
that longitudinal boosts and transverse rotations do not 
commute. Only in the rest frame |P| —> 0 do the relations 
for the longitudinal and transverse components coincide, 
owing to rotational symmetry. Interestingly, Burkardt [30j 
showed that the Ji relation (16) has a simple interpretation 
when the transverse component of the total AM is consid¬ 
ered in the rest frame. In order to obtain a momentum- 
independent relation for the transverse polarization, Ji, 
Xiong and Yuan [89] proposed to use the transverse com¬ 
ponent of the Pauli-Lubanski vector instead of the to¬ 
tal AM, but even in this case the separate quark and 
gluon contributions remain frame dependent [95119511951 
I63| . From now on, we will implicitly choose the z axis 


5 The direction of propagation of a laser beam provides a 
natural direction along which one can consistently define and 
measure separately the photon spin and OAM projections. 

6 This relation holds for a spin-1/2 target. For a spin-1 tar¬ 
get, a similar relation has been obtained by Taneja et al [US]. 


along the proton momentum and restrict ourselves to the 
longitudinal component of AM. 

The quark/gluon kinetic OAM is then obtained by 
subtracting the quark/gluon spin contribution, which is 
given in the MS scheme by the first Mellin moment of 
the quark/gluon helicity distribution 


( S q ) = \ J d xAq{x), 
( S G ) = J d xAG(x), 


(17) 


from the quark/gluon total kinetic AM 

(^ln) = (4 ? (n G )-(^ >G )- (18) 

Penttinen, Polyakov, Shuvaev and Strikman (PPSS) [140] 
showed that the quark kinetic OAM can also be directly 
expressed in terms of a two-parton twist-3 GPD 18711601 

[SUES] 

(LU = ~ JdxxG q 2 (x,0,0) (19) 

which enters the description of the longitudinal target spin 
asymmetry of deeply virtual Compton scattering |38| . The 
quark canonical OAM can also in principle be accessed via 
twist-3 GPDs by extracting the potential 0AM(Lp O t) and 
adding it to the kinetic OAM (L ^ in ), see Eq. (j9j. Explicit 
expressions can be found in Refs. mm- 


Since GPDs were shown to give access to AM and 
many transverse-momentum distributions (TMDs) would 
be identically zero in absence of parton OAM, there is 
some hope that TMDs could also give quantitative infor¬ 
mation about OAM. Based on quark model calculations, 
the following relations have been proposed [1251114311121 
[45] 

( L Ln) = [ dxd 2 k± [h\{x,k\) - gl L {x,k\)\ , 

, 2 ( 20 ) 

= - / dxd 2 fcj_ ^ h^(x,k\). 


Unfortunately, these relations cannot be derived directly 
from QCD and turn out to be valid only in a restricted 
class of models jl21U1191lll8| . In general, no direct quan- 
titave relations between OAM and TMDs should actually 
be expected as the former represents a correlation between 
parton position and momentum whereas the latter only 
provide information about the momentum distribution. 
Nevertheless, TMDs do provide some indirect information 
about the OAM content of the nucleon by constraining the 
wave function of the latter. 

Another possibility has also been explored. Burkardt m 
[29] suggested that the quark Sivers TMD f^(x, k\) and 
the quark GPD E q (x,t;,t) could be related by a chromo¬ 
dynamic lensing mechanism 


d 2 fc_L ^2/ir'Oufcl) 


oc / d 2 b ± i(b ± ) (S T x d b± ) z £ q (x,bl), 


(21) 
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where X(6j_) is called the lensing function and £ q (x, b\) = 
/ ( 2^)2 A 9 (:c, 0, — A\). This relation is supported 

by some model calculations [251134111311150] but can hardly 
be put on firmer theoretical grounds. A variation of it 
has however been used by Bacchetta and Radici pT3j to 
fit SIDIS data for the Sivers effect with the integral con¬ 
strained by the anomalous magnetic moments which led 
to a new estimate of (J£ in ) in good agreement with most 
common GPD extractions 


As shown by Lorce and Pasquini |1 1711122] , the most 
intuitive expression for the quark/gluon OAM is as a phase- 
space integral 



J dxd 2 k± d 2 b± (bj_ x k±) z p+%{x, fc_L, b±; W), 


( 22 ) 


where the Wilson line W ensures gauge invariance and 
the relativistic phase-space (or Wigner) distribution 
can be interpreted semi-classically as giving the quasi- 
probabilitjj^] for finding at the transverse position b± a 
quark/gluon with momentum (xP + , fcj_) inside a longitu¬ 
dinally polarized ( A' = A = +) proton. Unlike the version 
proposed in Refs. calm, this phase-space representation 
does not require any relativistic correction because it is 
based on the Euclidean subgroup of the light-front for¬ 
malism |144lf20t[50] . 

Under a Fourier transform, the phase-space distribu¬ 
tion is related to the generalized transverse-momentum 
dependent distributions (GTMDs) |132U121lll2()| . leading 
to the simple relation [ITHI661I85] 

(iw) = -J dxd 2 k ± ^F q f(x,0 1 k ± ,0 ± ;W). (23) 


The shape of the Wilson line determines the type of OAM 
[MmM]. For a staple-like Wilson line W^, like e.g. 
the one involved in the description of semi-inclusive deep- 
inelastic scattering and Drell-Yan process [37], Eq. (231 
leads to the canonical version of the quark/gluon OAM 
(-^can) = (Lyy ) irrespective of whether the staple is fu¬ 
ture or past-pointing [HHllH71IHTlllllllllfi| . For a straight 
Wilson line VV|, it leads to the kinetic version of the quark 
OAM (L^ in ) = (L q v ) [75] and the Ji-Xiong-Yuan (JXY) [ST 


definition of the gauge-invariant gluon OAM (Lf XY ) = 
(Ly i>l ) ■ Note that while the gluon kinetic OAM is given by 
Eq. ( |TT| ), the gluon JXY OAM is defined as 


L?x y = -J d 3 r E a i(r x V ab )A b J hys 

= tG , jG 

^can ' -^pot 5 

7 The relativistic phase-space distribution is generally not 
positive-definite because of the Heisenberg’s uncertainty prin¬ 
ciple. At best it can be seen as a weighing function in phase- 
space. Alternatively, one can consider the positive-definite 
Husimi distribution but one loses the probabilistic interpreta¬ 
tion of the corresponding position and momentum-space dis¬ 
tributions [57j . 


where 



7^ i pot . 


The JXY decomposition of the proton spin then reads 

J = S« + S G + LS XY + Lf XY -Ll at -L% t (26) 

with Ljxy = -^kin an d L lot = ~Lpot- 

While GTMDs provided for the first time a clear re¬ 
lation between canonical OAM and parton distributions, 
it is not known so far how to extract them from actual 
experiments, except possibly at small x m • Moreover, 
the GTMD F 14 does not reduce to any GPD or TMD, and 
so cannot be directly constrained. However, recent theo¬ 
retical developments [831184] and encouraging results [146] 
in the framework of large-momentum effective field theory 
(LaMET) [77] open the interesting possibility of comput¬ 
ing GTMDs and OAM directly on the lattice, see sec¬ 
tion |5.2| Moreover, effective models, constrained by ex¬ 
perimental data, also provide indirect access to GTMDs 
and OAM nT7irmHRHl[TT41ll3blll()0l . 


4 Current theoretical questions and puzzles 


From a theoretical point of view, it is naturally very inter¬ 
esting to check within models the different relations pre¬ 
sented in the previous section together with the AM sum 
rule. Several results have already been obtained in rel¬ 
ativistic quark models and perturbative models, but the 
situation is rather confusing and definitely requires further 
investigations. 

In Ref. m , Lorce and Pasquini compared, within two 
different relativistic light-front quark models, the value for 
the kinetic OAM obtained from Eq. (181 and the Ji rela¬ 


tion (161 with the value for the canonical OAM obtained 
from the model-dependent TMD relations (20) and the 


model-independent GTMD relation (23), see Table [lj All 
these three expressions agree on the total OAM, showing 
that the conservation of total AM is correctly implemented 
in the models. They however differ in the split of the lat¬ 
ter into the various flavor contributions. For the canonical 


OAM obtained from Eqs. (20) and (23), this can be at¬ 
tributed to the difference in the definition of the “center” 
of the target about which the OAM is computed [1 ID] . The 
difference in the flavor split between canonical and kinetic 
OAM is particularly puzzling because it is usually believed 
that the two forms of OAM should coincide in quark mod¬ 
els because of the absence of explicit gluon degrees of free¬ 
dom. Interestingly, a similar observation has been made 
earlier by Wakamatsu and Tsujimoto [ 159] within the chi¬ 
ral quark-soliton model based on the large-_/V c picture of 
QCD, and recently by Liu and Ma m within the axial- 
vector diquark model. So far no satisfactory explanation 
of these results has been proposed. 

In the scalar diquark model, the situation is also quite 
puzzling since, on the one hand, Burkardt and BC [33] 
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Table 1. Kinetic and canonical quark OAM obtained in the 
light-front constituent quark model (LFCQM) and the light- 
front chiral-quark soliton model (LFyQSM) [ 117] . 


LFCQM 

U 

d 

Total 

(^kil^Eq. p] 

0.071 

0.055 

0.126 

{-^can^Eq. j20| 

0.131 

-0.005 

0.126 

(-^can^Eq. \23) 

0.169 

-0.042 

0.126 

LFxQSM 

u 

d 

Total 

( L kin)Eq. {Ts| 

-0.008 

0.077 

0.069 

(-^- / can)Eq. j20| 

0.073 

-0.004 

0.069 

(-^can)Eq. \23) 

0.093 

-0.023 

0.069 


sity is converted into an additional contribution to OAM 
density plus a superpotential term. One has to be careful 
when defining AM at the density level as superpotentials 
modify the distribution and in turn destroy the interpre¬ 
tation as AM density. It is therefore crucial to keep track 
of these superpotential terms when discussing AM at the 
density level. 

A first approach has been proposed by Hagler and 
Schafer [58] where the higher moments of the canonical 
AM densities in the light-front gauge A + = 0 are defined 
by a tower of generalized canonical AM operators, like e.g. 

^can 1 = / d 3 r 4>'Y + ( r ± x id±) z iD^ 1 ■ ■ ■ (27) 


found that the quark canonical and kinetic OAM do co¬ 
incide using the Pauli-Villars regularization while, on the 
other hand, Liu and Ma [ 107] found that they differ using 
a form factor regularization. For an electron in QED to 
order a , the difference between electron canonical and ki¬ 
netic OAM is non-zero, as expected for a gauge theory. It 
turns out to be finite but the precise value depends on the 
regularization prescription 


108] . Since this difference 
is attributed to the potential OAM ([9|, Liu and Ma 1108] 
computed explicitly the latter and found the suprising re¬ 
sult that it does not account for the observed difference. 
Lowdon m emphasized that surface terms usually as¬ 
sumed to vanish in field theory could actually contribute. 
Liu and Ma then computed explicitly the contribution 
from a surface term and found it to be non-zerc]f] Unfor¬ 
tunately, the combination of potential OAM and surface 
term is still not able to explain the difference between 
canonical and kinetic OAM, irrespective of the adopted 
regularization prescription. 

Obviously, further investigations are needed to clar¬ 
ify the situatio n. In particular, explicit calcu lati ons of the 
PPSS relation (191 and the GTMD relation (231 with dif¬ 
ferent Wilson lines should shed some light. Also, special 
attention to Lorentz covariance, equations of motion, reg¬ 
ularization prescription and surface terms should help re¬ 
solve the current puzzles. 


So far, we have focused our discussions on the AM in¬ 
tegrated over all phase space. It is however possible to 
define densities of AM, and therefore map the latter in 
position and momentum space. There is however a lot 
of confusion due to the large number of definitions pro¬ 
posed in the literature. One of the reasons for this confu¬ 
sion is that densities differing by a so-called superpotential 
term d a X^ a ^ '(r), where square brackets mean antisym- 
metrization of indices, lead to the same integrated quanti¬ 
ties. The archetypical example is the total AM of a system, 
which can be obtained either from the canonical density 
or the Belinfante-improved density. In the former case, the 
AM density is naturally decomposed into an OAM density 
and a spin density, while in the latter case, the spin den- 

8 This is actually expected since plane waves are somewhat 
pathological as they do not vanish at infinity, see [941 . 


Hoodbhoy, Ji and Lu [69] adopted the same strategy to 
define the corresponding higher moments of the kinetic 
AM densities and concluded that 


(L q kin )(x) = f lW(x,0,0) + E«(x,0,0)} - \Aq{x). (28) 

Note however that (L%. in )(x) is not defined from the simple 
tower 

fq,+—+ __ 

^kin — 

J d 3 r'i/jj+(iD + y(r± x iD±_) z (iD + ) n ~ 3 ^ ^ 

3=0 J 


but involves another complicated tower AlA *" such that 
the sum L^Sn + = "' + + AL + + evolves as a leading- 

twist operator. Moreover, as noted in Ref. J5S], Hoodbhoy, 
Ji and Lu assumed that the surface terms vanish for all 
the Mellin moments, and therefore implicitly assumed that 
the superpotential relating Belinfante and kinetic AM op¬ 
erators does not affect the x-dependence. In conclusion, 
the interpretation of the RHS of Eq. (28) as the density 
of quark kinetic OAM should be considered as an ad hoc 
definition rather than following from a natural operator 
definition. 

The alternative derivation of the Ji relation for quarks 
in the rest frame proposed by Burkardt [30] suggested that 
the integrand | [H q (x, 0,0) + E q (x, 0, 0)] could naturally 
be identified with an angular momentum decomposition of 
the transverse component of the quark angular momentum 
in a transversely polarized target. Ji, Xiong and Yuan tried 
to justify further this partonic interpretation by consider¬ 
ing a slight variation of Burkardt’s derivation E23 and by 
extending it to gluons [5TJ. However, a careful inspection 
reveals that this interpretation is not well founded. Indeed, 
the simple partonic interpretation proposed by the authors 
is based on the “leading-twist” part f d 2 r± r±T ++ of the 
Lorentz transformation generator f d 2 r± M ++± . Lorentz 
invariance is then invoked to include the “higher-twist” 
part f d 2 r±r + T +± . While this argument is perfectly valid 
for the x-integrated quantities, it manifestly cannot be 
applied at the density level. Moreover, as pointed out in 
Refs. [2152], Ji, Xiong and Yuan incorrectly interpreted 
M ++± as the transverse angular momentum operator. In 
light-front quantization, M ++± corresponds actually to 
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the transverse boost operator which is kinematical and 
therefore has a simple partonic interpretation, unlike the 
transverse rotation operator M -1 ^ which is dynamical 
and hence complicated due to the presence of interac¬ 
tion terms [25] . In order to obtain a boost-invariant sum 
rule, Ji, Xiong and Yuan [ST] focused on the transverse 
components of the Pauli-Lubanski vector and therefore 
included the higher-twist contribution M H but the 
strategy remained the same: do away with the higher- 
twist contribution using Lorentz symmetry [52]. More¬ 
over, the contribution associated with the C q,G (0) energy- 
momentum form factor has been arbitrarily discarded by 
Ji, Xiong and Yuan mm, based on the argument that it 
is higher-twist and does not contribute to the total Pauli- 
Lubanski vector. This contribution is however important 
once one considers quarks and gluons separately mm 
1631195] . Moreover, energy-momentum form factors being 
Lorentz scalars, naive twist counting of amplitudes does 
not provide any information about their magnitudes and 
therefore does not justify neglecting C q ’ G ( 0) in front of 
the other energy-momentum form factors. 

Clearly, the interpretation of the integrand of the Ji 
relation (16) as the density of parton AM is not justified. 


(■ LlJ(x)=x[H q (x,0,0)+E q (x,0,0 ) 

+G\ 0,0,0) - 2 G\(x, 0,0)] - Aq(x) 


(30) 


Harindranath and Kundu [61] defined a natural density of 
quark and gluon canonical OAM in the light-front gauge. 
Bashinsky and Jaffe m then improved this definition by 
providing an expression invariant under residual gauge 
transformations. The density of parton OAM can in fact 
be defined in phase space and is easily expressed in terms 
of gauge-invariant phase-space distributions HT711122Hlll| 

{L^)(x, k±,b ± ) = (b_ l x kj_) z p q ^ G (x,k±,b ± ;W). (32) 


Like in the integrated version (22), a staple-like Wilson 


This is actually expected since the presence of the covari¬ 
ant derivative in the kinetic definition of OAM necessarily 
implies the contribution of higher-twist parton distribu¬ 
tions. Hagler, Mukherjee and Schafer ESI proposed the fol¬ 
lowing density of OAM involving both twist-2 and twist-3 
GPDs 


which when combined with Eq. (28) led them to another 
form 

(L q ki „)(*) = -x [G q 2 (x, 0 , 0 ) - 2 G\(x, 0 , 0 )]. (31) 

The interpretation of Eq. ( filOj ) as the distribution of OAM 
has however been obtained within the Wandzura-Wilczek 
approximation, where the distinction between canonical 
and kinetic OAM disappears. Ji, Xiong and Yuan [81] de¬ 
fined the densities of quark/gluon kinetic and potential 
OAM via moments of the form (29) and showed how the 
latter are related to D and F-type twist-3 GPDs. One 
has however to keep in mind that the gluon contributions 
defined in this way are those entering the JXY decompo¬ 
sition (26) and not the standard kinetic decomposition. A 


line leads to the canonical version of the quark/gluon 
OAM, but this time it depends on whether the staple 
points toward the future or the past W c . Indeed, 
based on the constraints imposed by parity and time- 
reversal symmetries, the relativistic phase-space distribu¬ 
tion of unpolarized partons inside a longitudinally polar¬ 
ized nucleon can be decomposed into four modes 

P Q ++ = P\' G + (b± ■ k ± ) p\ G + (b ± x k ± ) z pt 
+ (b ± xk ± ) z (b ± -k ± )pf G , 

where pf G = pf G (x, k']_, (bj_-fcj_) 2 , by, W). The functions 
pi and p 2 are associated with the real and imaginary parts 
of the GTMD Fn, whereas p% and p 4 are associated with 
the real and imaginary parts of the GTMD Fi 4 . The fac¬ 
tor (bj_ • k±) indicates that P 2 and p\ are naive T-odd, i.e. 
change sign under VW, 0 yV c . Note also that the only 
contribution to the canonical OAM surviving integration 
over b± or k± is pz- For a straight Wilson line W), one 
obtains the density of JXY quark and gluon OAM [791181] 
provided that one integrates at least over k± m- 

We close this section with a discussion of the Ji re¬ 


more detailed discussion of the connection between twist- 
3 GPDs and densities of AM has been given by Hatta and 
Yoshida m- They stressed in particular that the den¬ 
sity of kinetic OAM is ambiguous because D-type twist-3 
GPDs involve two fractions of longitudinal momentum x\ 
and X 2 ■ This ambiguity does not show up for the density 
of canonical OAM since the latter has support only at 
Xi = X 2 - This boils down to the fact that, contrary to or¬ 
dinary derivatives, covariant derivatives do not commute 
and therefore do not admit a unique non-local generaliza¬ 
tion mi- 

Motivated by the simple partonic picture behind the 
canonical AM decomposition of Jaffe and Manohar [72] . 


lation (16) in position space. Since the t dependence of 
twist-2 GPDs contains the information about the spatial 
distribution of partons [251128] . it has been suggested by 
Polyakov [1421151] that the Ji relation generalized to t 7 ^ 0 
contains the information about the spatial distribution of 
Belinfante AM 

( J kkf )(*)=! J dxx[H q ’ G (x,0,t) + E q ’ G (x,0,t)}. (34) 

Within the scalar diquark model with Pauli-Villars regu¬ 
larization where (F^ in ) = (L q an ), Adhikari and Burkardt [J 
observed that the quark kinetic and canonical spatial den¬ 
sities do not coincide (L kin )(b±) / (L q an )(b±). Note how¬ 
ever that these spatial densities have been defined in the 
Drell-Yan frame as simple two-dimensional Fourier trans¬ 
form (L q ){b ± ) = f < 0rk e ~ iA± b± {L q )(t = -A\). While 
this is justified for the canonical OAM and helicity distri¬ 
butions, the actual spatial density of total AM is given by 
the Fourier transfornj®] of {J kin ){t) + t ]|(J kin )(t) E23- The 

9 In the non-relativistic interpretation based on the large- 
N c picture used by Polyakov, the spatial distribution is given 
by the three-dimensional Fourier transform of (J^ in ){t) + 
T Il42l . Note that a quadrupole contribution has 

implicitly been discarded in Ref. [51] by doing the substitution 

-A 3 , 
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second term does not contribute to the total AM since 
the integral over all space amounts to set t = 0. Note also 
that the spatial density of quark kinetic OAM has been 
defined as (Ll in )(b±) = (</£„) (b±)-(S 9 ) (bjJ, but this re¬ 
lation is actually spoiled by the superpotential term relat¬ 
ing Belinfante and kinetic forms of the energy-momentum 
tensor [951 . Further investigations are once more needed. 

5 The orbital angular momentum from the 
lattice 


0 [152| . This EMT can be decomposed into quark and 
gluon kinetic contributions 

= TZ, q + TZ, G ’ ( 37 ) 

where and G have the following form 

rz, q =! E WB" + YD^ f , (ss) 

/ 

T&l\ G = -2 Tr[F^F v x ] + \ g^ Tr[F 2 ]. (39) 


As mentioned in the Introduction, the quark spin as deter¬ 
mined from DIS experiments is about 25% of the proton 
spin and recent global analyses [TUI 139| including high- 
statistics 2009 STAR [T] and PHENIX [5] data show an 
evidence of non-zero gluon helicity in the proton. For Q 2 = 
10 GeV 2 , it is found that the gluon helicity distribution 
AG(x, Q 2 ) is positive and different from zero in the mo¬ 
mentum fraction range 0.05 < x < 0.2 and the integral 
J 0 ' 05 dx AG(x) is 0.17 ± 0.06 - about 35% of the proton 
spin |l39j . However, the results have very large uncertainty 
in the small x-region [mum Since the sum of the quark 
and gluon spin do not saturate the total proton spin, it 
leaves quite a bit of room for the OAM of quarks and 
gluons. 

Now that we have discussed the decomposition of the 
proton spin in terms of the quark and gluon contributions 
in Eqs. (10) and (13), and how to access OAM using GPDs 
and GTMDs, it would be useful and imperative to calcu¬ 
late them within lattice QCD which consists in ab initio 
calculations with controllable statistical and systematic 
errors. 


5.1 Mellin moment approach 


We shall first consider Ji’s decomposition (131 which in¬ 


volves the quark spin, quark kinetic OAM, and gluon to¬ 
tal AM. Since the decomposition of longitudinal AM is 
boost-invariant, one can calculate the Mellin moments of 
the GPDs H and E through the matrix elements of the 
energy-momentun tensor (EMT) near the rest frame like 
in the case of the parton structure functions. The lattice 
calculation of physical quantities involving small lattice 
momenta near the lowest non-zero lattice momentum, i. e. 
pa ~ ITT with a the lattice spacing and La the physical 
lattice size, are less noisy than those involving pa ^> jT . 

The Lorentz group generator J is given by [162] 


J ^ = J d 3 rM°^ v (r). 


(35) 


Here M 0 ^ is the Lorentz group generator density which 
is defined as 


M^ a (r) = T^r" - T^r 01 , (36) 


with D v =D V — D". The AM operator is defined from the 
Lorentz group generator in Eq. (351 as 


J 1 = 


1 ^ijk 


d 3 r M°i k (r). 


(40) 


Since the quark kinetic OAM L^ in in Eq. ([2]) and the 
gluon total AM operators Jjf in in Eq. (14) depend on the 


spatial vector r, a straightforward lattice calculation of the 
operators with forward matrix elements are complicated 
by the periodic condition of the lattice, and may lead to 
wrong results [ 163] . Instead of calculating L^ in and J^ in 
directly, they can be obtained from the energy-momentum 
form factors in the nucleon which parametrize off-forward 
matrix elements of the EMT [mimiTTB] . 


To carry out lattice calculations, one works with the 
Monte-Carlo ensembles of the 2- and 3-point path-integral 
correlation functions in the Euclidean space. The Euclidean 
energy-momentum operators for quarks and gluons are [40] 

T 4 t q = + Tz^Bf]?/’/, (41) 

/ 

Tla = 2iTr[F E a F E \, (42) 

where the Euclidean 7 ^ matrices satisfy { 7 ^, 7 ^} = 2 5^ 
and the Euclidean components of a four-vector are defined 
as = ( p,ip °). In the following, the Pauli-Sakurai rep¬ 
resentation for the gamma matrices |1331fTBl is used and 
the covariant derivative is the point-split lattice operator 
involving the gauge link [89]. 

The form factors of the quark and gluon EMTs are 
defined as 


ip', s'\T£ a \p, S) = 4 u E (p\ s’)rl a u E {p , s) 

with a = q,G ancp*] 

AE = (A4^+7^4)T“(-A 2 ) 

{Pi&ia T Pi(J 4a )A a a 2 


(43) 


2 m 
T 3 “(-Z\ 2 ), 


T“(-A 2 ) (44) 


iA 4 Ai 2 


where is the Belinfante-improved EMT which is gauge 
invariant, symmetric = T M " and conserved = 


10 Note that the energy-momentum form factor T 4 does not 
contribute to (44 1 as it is associated with the structure An = 0. 
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where P = p f p , A = p' — p, and the normalization con¬ 
ditions for the nucleon spinors are 

u E (p , S) u E (p , S) = 1 

. it) + M 

Y^, u {P,S)u E {p,S)= 2M . (45) 

The EMT form factors T\ and T 2 resemble the Dirac and 
Pauli form factors Fi and F 2 associated with the electro¬ 
magnetic current. While the nucleon electric charge and 
magnetic dipole moment respectively equal F\ and F 1 +F 2 
at A 2 = 0, the quark and gluon kinetic momenta and total 
kinetic AM fraction are shown to be respectively equal to 
Tf and Tf + Tf at Z \ 2 = 0 [74|. 


Lattice theory has been able to accommodate chiral 
symmetry of the vector gauge theories finally, the lack of 
which has hampered the development of chiral fermions on 
the lattice for many years. It is shown that when the lat¬ 
tice massless Dirac operator satisfies the Gingparg-Wilson 
relation 'y E D E + D Er yf = aD E "/ E D E with the overlap 
fermion being an explicit example |138j , the modified chi¬ 
ral transformation leaves the action invariant and gives 

rise to a chiral Jacobian factor J = e _2lQTr [75 { - 1 ~ 2 aD 
from the fermion determinant jl24i . The index theorem [64j 
shows that this Jacobian factor carries the correct chiral 
anomaly. Thus, the Atiya-Singer theorem is satisfied on 
the lattice with a smooth gauge configuration, i. e. the to¬ 
tal topological charge equals the index (number of left- 
handed zero modes minus right-handed ones) of the con¬ 
figuration. Further, it is shown that the local version of 
the overlap Dirac operator gives the topological charge 
density operator in the continuum HMUmZ!, namely 

Tr[7jf (1 — \ aD E v (r,r))\ = g^T \[F E v F E ^) + 0{a). (46) 

This definition is superior to previous lattice constructs 
with gauge links. Since the overlap operator is exponen¬ 
tially local, the topological charge operator so defined is 
chirally smoothed with neighboring links. This operator 
has been used to study the topological properties of the 
vacuum, and charged coherent membranes with lower¬ 
dimensional long-range order mm were discovered which 
have implications on chiral symmetry breaking. 

In view of the fact that the overlap operator definition 
of the topological charge has good properties, a similar 
connection between F' E v and the overlap Dirac operator 
has been derived mm 


Tr. 


[<V- D ov( r >r)] = c T a F^ 




(r) + 0(a 3 ), (47) 


where Tr s is the trace over spin and ct = 0.11157 is the 
proportionality constant in the continuum limit for the 
parameter k = 0.19 in the Wilson kernel of the over¬ 
lap operator which is used in the lattice calculation (1041 . 
The overlap Dirac operator D E v (r,r') is exponentially lo¬ 


cal and the gauge held F E V as defined in Eq. ( |47| is chi¬ 


rally smoothed so that it admits good signals forthe gluon 
momentum and total AM in the lattice calculation I4T1. 


This method proved to be much less noisy than the stan- 
drad approach based on the clover-leaf definition from the 
gauge links m- 


5.1.1 Sum rules and renormalization 


The momentum and total AM contributions from quarks 
and gluons depend on the renormalization scale and scheme 
individually, but their sums do not because the total mo¬ 
mentum and AM of the nucleon are conserved. One can 
use these sum rules as the renormalization conditions on 
the lattice. 

Substituting the EMT matrix elements in Eq. (431 to 
the matrix elements which define the AM in Eq. (401 and 
similarly for the momentum, it is shown m that 


(4„ G > = Z e g I [Tf G (0) + Tf G (0)], (48) 

(x) q , G = Z E G T?' G ( 0), (49) 


where Z qG is the renormalization constant for the lattice 
quark/gluon operators, and (x) q> G is the second Mellin 
moment of the unpolarized parton distribution function 
(PDF) f dxx fi’ G (x), which represents the (average) mo¬ 
mentum fraction carried by quarks or gluons inside a nu¬ 
cleon. The other form factor, Tf G (0), can be interpreted 
as the anomalous gravitomagnetic moment in analogy to 
the anomalous magnetic moment F 2 ( 0) [1501122) . 

Since momentum is always conserved and the nucleon 
has a total spin of 1/2, we write the momentum and AM 
sum rules using Eqs. ([48| and (49), as 


Z e T?(0) + Z e T g (0) = 1, (50) 
{Z L q [T/(0) + T|(0)] + Z E [T g (0) + T g (0)]} = (51) 


It is interesting to note that from Eqs. (50) and (511, one 
obtains that the sum of the T 2 (0)’s forthe quarks and 
gluons is zero, i.e. 


ZfT|(0) + ZfTf(0) = 0. 


(52) 


The vanishing of the total anomalous gravitomagnetic mo¬ 
ment (T 2 (0) = 0) in the context of a spin-1/2 particle was 
first derived classically from the post-Newtonian manifes¬ 
tation of equivalence principle ( 88 ]. More recently, this has 
been proven by Brodsky et al. |22] for composite systems 
from the light-front Fock space representation. 

One can use these sum rules and the raw lattice results 
to obtain the lattice renormalization constants Z E and 
Z G , and then use perturbation theory m to calculate 
the quark-gluon mixing and renormalization in order to 
match to the MS scheme at 2 GeV which preserves the 
sum rules. 


5.1.2 Results of lattice calculations 

Before presenting the lattice results, we should point 
out that the 3-point functions for quarks which are needed 
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(a) 


(b) 


Fig. 1. Quark line diagrams of the 3-point function with 
current insertion in the Euclidean path integral formalism, 
(a) Connected insertions (Cl), and (b) disconnected insertions 
(DI). 


to extract the form factors in Eq. (44) have two topologi¬ 
cally distinct contributions in the path-integral diagrams: 
one from connected insertions (Cl) and the other from 

Fig-0 


disconnected insertions (DI) |H)lUlU^U99Ull)bj . see 
They arise in different Wick contractions, and it needs to 
be stressed that they are not Feynman diagrams in per¬ 
turbation theory. In the case of Cl, quark/antiquark fields 
from the operator are contracted with the quark/antiquark 
fields of the proton interpolating fields. It represents the 
valence and the higher Fock space contributions from the 
Z-graphs. In the case of DI, the quark/antiquark fields 
from the operator contract themselves to form a current 
loop, which represents the vacuum polarization of the dis¬ 
connected sea quarks. Also, although the quarks lines in 
the loop and the nucleon propagator appear to be dis¬ 
connected in Fig0)b), they are in fact correlated through 
the gauge background fluctuations. In practice, the un¬ 
correlated part of the loop and the proton propagator is 
subtracted. For the nucleon, the up and down quarks con¬ 
tribute to both Cl and DI, while the strange and charm 
quarks contribute to the DI only. 

As for the quark OAM, lattice calculations have been 
carried out for the Cl pnilMlIMlIMlETlITinimUTT] . They 


are obtained by subtracting the quark sp in c ontributions 
from those of the quark kinetic total AM ( fl8| ) . It has been 
shown that the contributions from u and d quarks almost 
cancel each other. Thus for Cl, quark OAM turns out to 
be small in the quenched calculations |130ll551 and nearly 
zero in dynamical fermion calculations p4ll2lHT0lll481lllj . 

The first complete decomposition of proton spin into 
quark spin, quark OAM and gluon total AM together with 
quark and gluon momenta has been carried out by Deka 
et al. (\QCD Collaboration) [2T]. The calculation is done 
on a quenched lattice with 3 valence quark masses and 
extrapolated to the physical pion mass. The results are 
presented in Tableland Fig. [2] 

We see from Table [2] that the strange momentum frac¬ 
tion (x) s = 0.024(6) is in the range of uncertainty of 
(x) s from the CTEQ fitting of the PDF from experiments 
0.018 < ( x) s < 0.040 [HI]. The gluon momentum fraction 
of 0.334(55) is smaller than, say, the CTEQ4M fit of 0.42 
at fj, = 1.6 GeV [90), but only by 1.5cr. The smallness of 
this value of ( x)g in comparison to the experiment could 
be in part due to the fact that this is a quenched calcu- 



(a) 




(c) 


Fig. 2. Pie charts for the quark and gluon (a) momentum 
fractions, (b) AM fractions, and (c) spin and OAM fractions. 
The left panels show the quark contributions separately for 
Cl and DI, and the right panels show the quark contributions 
for each flavor with Cl and DI summed together for u and d 
quarks. 


lation. It is expected that the gluon momentum fraction 
will be larger than the present result when dynamical con¬ 
figurations with light fermions are used in the calculation. 

For the unrenormalized lattice results, [?2 i (0)-l-T^(0)]ci 
is found positive and T§ (0) negative, so that including the 
small [T^(0) +T^(0) +T|(0 )]di, the total sum can be nat¬ 
urally constrained to be zero, see Eq. (52) with the lattice 
normalization constants Z^ = 1.05 and Zq = 1.05 close to 
unity. As discussed earlier, the vanishing of the total T 2 (0) 
is a consequence of momentum and AM conservation. 

The flavor-singlet g° A = g\ + g d A + g s A = AS = 2(S' 9 ), 
which represents the quark spin fraction to the nucleon 
spin, has been calculated before on the same lattice pFfl . 
The quark spin contribution are subtracted in Table[2]and 
Fig. 0 from the quark total AM contributions (j£~) to 
obtain the quark OAM contributions (L^. in ). As we see in 
Table 0 the OAM fractions 2(L£ in ) for the u and d quarks 
in the Cl have different signs and add up to 0.01(10), 
i.e. essentially zero. This is the same pattern which has 
been seen with dynamical fermion configurations and light 
quarks, as pointed out earlier. The large OAM fractions 
2(Tkin) for the u/d and s quarks in the DI is due to the fact 
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Table 2. Renormalized results for connected (Cl) and disconnected (DI) insertions in MS scheme at g = 2 GeV. 



CI(-u) 

CI(d) 

CI(w + d) 

DI (u/d) 

DI(s) 

Gluons 

( X } 

0.413(38) 

0.150(19) 

0.565(43) 

0.038(7) 

0.024(6) 

0.334(55) 

t 2 ( 0 ) 

0.286(108) 

-0.220(77) 

0.062(21) 

-0.002(2) 

-0.001(3) 

-0.056(51) 

2 (Jkm) 

0.700(123) 

-0.069(79) 

0.628(49) 

0.036(7) 

0.023(7) 

0.278(75) 

9A 

0.91(11) 

-0.30(12) 

0.62(9) 

-0.12(1) 

-0.12(1) 

- 

2 (I/kin) 

-0.21(16) 

0.23(15) 

0 .01(10) 

0.16(1) 

0.14(1) 

- 


that g\ 4 in the DI is large and negative, about —0.12(1) 
for each of the three flavors. All together, the quark OAM 
constitutes a fraction of 0.47(13) of the nucleon spin. The 
majority of it comes from the DI. 

As far as the spin decomposition is concerned, it is 
found that the quark spin constitutes 25(12)% of the pro¬ 
ton spin, the gluon total AM takes 28(8)% and the rest is 
due to the quark kinetic OAM which is 47(13)%. 

Since this calculation is based on a quenched approxi¬ 
mation which is known to contain large uncontrolled sys¬ 
tematic errors, it is essential to repeat this calculation 
with dynamical fermions of light quarks and large physi¬ 
cal volume. However, we expect that the quark OAM frac¬ 
tion may still be large in the dynamical calculation. First 
of all, the quark spin fraction is about 20% - 30% from 
DIS and the gluon momentum fraction from CTEQ4M 
fit gives (x)g = 0.42 at g = 1.6 GeV (90j which is not 
far from our calculation at /i = 2 GeV. Since the gluon 
total AM is (x)g + Tp(0), if we assume that T| ? (0) in 
the dynamical calculation is not too far from the small 
Tjp(0) = —0.056(51) obtained in the quenched approxi¬ 
mation, the gluon total AM will be ~ 36%. In this case, 
the quark OAM fraction will still be at the 40-50% level. 
This will be checked by the upcoming lattice calculations 
with dynamical fermion configurations at physical pion 
mass. 

In the naive constituent quark model, the proton spin 
comes entirely from the quark spin. On the other hand, 
in the Skyrme model [5] the proton spin originates solely 
from the OAM of the collective rotational motion of the 
pion field Ez]. What is found in the present lattice calcu¬ 
lation seems to suggest that the QCD picture, aside from 
the gluon contribution, is somewhere in between these two 
models, indicating a large contribution of the quark OAM 
due to the meson cloud (qq pairs in the higher Fock space) 
in the nucleon. 


5.1.3 Asymmetric energy-momentum tensor 


Let us stress that the standard Mellin moment approach 
is based on the Belinfante-improved EMT which is sym¬ 
metric and therefore leads only to the total AM of quarks 
and gluons, see Eqs. (351 and (36). In the quark sector, 
one can alternatively use the asymmetric EMT 


= ( 53 ) 

/ 


which differs from the Belinfante-improved EMT (38) by 
a gauge-invariant antisymmetric superpotential term [95] 


TZ, g = rz, q -1 d a (e^ 3 Y, ( 54 ) 

/ 


provided that the QCD equations of motion = 

0 hold. From this asymmetric EMT, one can directly com¬ 
pute the quark kinetic OAM 


L 


kin,g 


— 1 ijk 
~ 2 t 


d 3 T- (T.°. k T-? 


nOj 


r l J kin,^" - T^ q r- 


(55) 


The parametrization of the off-forward matrix elements of 
the asymmetric EMT, beside the structures in Eq. 
involves a new structure ( 74 Pi — J 1 P 4 ) T|(—Zl 2 ) [14111161 . 
The quark spin and kinetic OAM fractions are then re¬ 
spectively given by — T| and T q + T| + T| at A 2 = 0 [131 
mm- To the best of our knowledge, this approach has 
not been investigated so far on the lattice. In particular, it 
would be interesting to compare the lattice calculation of 
—and g q A since both represent the quark spin fraction 
at A 2 = 0 obtained from two different operators. 


5.2 Large momentum approach 

To carry out lattice calculations of all the terms in the 
generic sum rule of Eq. ([!]) , one needs to define the fields 
A£ ure and A^ hys on the lattice. Moreover, the partonic 
picture of OAM is naturally depicted in the light-front 
formalism and related to GPDs and GTMDs as discussed 
in Sec. [3j Unfortunately, the light-front coordinates are 
not accessible to lattice QCD calculation since the latter 
is based on Euclidean path-integral formulation. 

To bridge the gap between the light-front formulation 
and the lattice calculation, it is suggested by Ji m that 
light-front parton distributions can be accessed as the in¬ 
finite momentum limit of equal-time correlators. To do 
so, one calculates the equal-time correlators on the lattice 
at large but finite hadron momentum, called quasi parton 
distributions, and match them to the MS scheme at a cer¬ 
tain scale g in the continuum. After the usual continuum 
extrapolation of the lattice results, a large momentum ef¬ 
fective field theory (LaMET) [78], which takes care of the 
non-commuting UV and P z —>• 00 limits, is used to match 
these quasi parton distributions to ordinary parton distri¬ 
butions on the light-front [771ll66UT2i61 . 

This approach can be adopted to many observables, 
but one needs to first verify that the operator used in 
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the lattice calculation, when boosted to the infinite mo¬ 
mentum frame (IMF), coincides with the one on the light- 
front. The gluon lrelicity distribution is defined by Manohar 
in terms of the light-front correlation function 1129 


AG{x)S+ = 


2 xP+ 


dz~ 
2i r 
iab 


tP + z~ 


(P, S\F +a ’ a (0)W ab {0, z~)F + a ' b {z~)\P , S), 


(56) 


where W(Q,z~) = V[e~ 19 fo d C A+ C -° x )] is a light-like 
Wilson line. After integrating over the longitudinal mo¬ 
mentum fraction x, this light-front correlator reduces to 
the proper definition of gluon light-front helicity |65 l l8.‘l l . 
Furthermore, it has recently been shown by Ji, Zhang and 
Zhao [331 that this quantity coincides with the gluon spin 
contribution defined by Chen et al. [35][36l in the IMF with 
Aphys satisfying the non-Abelian Coulomb condition [3(3] 

V ■ A phys = 0. (57) 


gauge and the Coulomb gauge-fixed potential A c (r) it¬ 
self. Therefore, one can construct the quark and gluon 
canonical OAM L q an and n as local operators on the 
lattice from U£ ure (r) and Ap hys (r). This strategy has re¬ 
cently been applied to calculate the gluon light-front he¬ 
licity ag manso]. 

Note that L q an resembles the kinetic OAM in Eq. ([2]), 
except that the covariant derivative D M in the latter is 
replaced by D£ ure . By analogy with Eq. ( |38| ), one might 
think to get L q an by subtracting the quark spin from the 
quark total AM obtained from the symmetric canonical 
EMT 




+ 7 1 ibf 

pure 1 / -^pureJrv 


(62) 


However, unlike the kinetic AM case in Eq. (36), the den¬ 


sity of total AM is not simply related to this symmetric 
canonical EMT 


In other words, the longitudinal gluon spin operator turns 
into the helicity operator in the IMF. Similarly, Zhao, 
Liu and Yang m have shown that the quark and gluon 
canonical OAM defined by Chen et al. [351136] with the 
condition (57) in IMF coincide with the expressions de¬ 
rived explicitly by Hatta [6H] and Lorce inn with staple 
light-like Wilson lines, namely 


lim (L 

P z —>oo 


q,G\ 

can / 


Aphys=0 



(58) 


To carry out a lattice calculation of the canonical OAM 
and gluon spin, one possibility is to define A(j ure and Ap hys 


with the constraint (57) on the lattice. It is shown in 


Refs. [16911167] that their solutions are related to the gauge 
link (7 M (r) which connects the site r to r + ag (a is the 
lattice spacing). After the gauge link is fixed to that of 
the Coulomb gauge on the lattice with a gauge rotation 
matrix g c 


U^(r) = g c {r)U£{r)g c 1 (r + ajl), 


one can define the links U£ uie and t/p hys as 


(59) 


^pureO) = gc{r)g c (r + a£), 

Khy^) = - Kure(r), 


(60) 


so that 


K^) = i- g [U^ re {T)-t], 

^ph ys ( r ) = ffc(r)A^(r)g- 1 (r) + 0{a). 

It is straightforward to prove [16711169111 TO] that not only 
A£ ure and Ap hys so constructed satisfy the gauge transfor¬ 
mation properties in Eqs. © and (pj), but the pure gauge 
condition in Eq. Q and the non-Abelian Coulomb con¬ 
dition in Eq. ( [57[ are also satisfied to G(a). Note that 
A£ ure and A 9 hys are solutions which are obtained from 
the gauge rotation matrix g c (r ) which fixes the Coulomb 


M^ a ( r) 4 7^° r v - r“, 

can,gv' ) / 7 can,<7' 7 can,<7' j 


(63) 


but involves extra terms. Therefore, there is no such a 
relation as Eq. (48) in the canonical case. If, however, one 
adopts the asymmetric canonical EMT, 


/ 


(64) 


one does obtain, similarly to the kinetic case in Eq. (55), 
directly the quark canonical OAM 


jJ — 1 e ijk 


^ 3 ^ (^rOk 


r (T ok r 9 — r k 1 

' V^can ,q ^can,o' /' 


(65) 


Because of the Stueckelberg/background dependence of 
the decomposition of the gauge potential <[3j) , more struc¬ 
tures in the parametrization of the off-forward matrix el¬ 
ements of T££ nq are allowed. In particular, for the natu¬ 
ral choice determined by the QCD factorization theorems, 
one has to introduce a light-like vector nF. Once general¬ 
ized form factors are separated in lattice QCD calculations 
through proper linear combinations of matrix elements at 
different kinematics, then the quark canonical OAM can 
be obtained through the relation [1T6] 

(L q am )=At’ q ( 0,0). (66) 

Other similar relations have also been obtained for the 
gluon canonical OAM and the quark/gluon spin contribu¬ 
tions. 

Alternatively, one can directly calculate the canonical 
OAM from the off-forward matrix element, just like the 
magnetic moment can be obtained from the off-forward 
spin-flip matrix element [3] , 


e ljk (p,S\ J d 3 r , tp'y 0 r J iD k uie 'ip\p, S) = 

e i9k dm (P',5'1 J d 3 r^°^e^ r D k UI ^\P,S). 


(67) 
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The challenge here is to have a reliable extrapolation to 
the A — > 0 limit which can be ameliorated with a large lat¬ 
tice or twisted boundary condition to have data at small 

A. 


The third option is the calculation of quark and gluon 
OAM from GTMDs. Exploratory lattice calculation of 
TMDs have been carried out in Refs. (13611137(1461147] where 
a staple-like Wilson line was used in the calculation 
and the result were evolved to the IMF to obtain the 
light-front correlator. This approach can be generalized to 
GTMDs and use Eq. (231 to obtain the quark and gluon 
canonical OAM. Since the link between the bilinear quark 
fields is non-local, one of the challenges in this approach 
is the renormalization of the non-local operator on the 
lattice. 


6 Effective theory of baryons 


Many estimates of quark spin and OAM contributions of 
the nucleon are based on quark models. However, quark 
models are not realistic effective theories of QCD, since 
they do not have chiral symmetry, a salient feature of 
QCD whose dynamics governs light-quark hadron struc¬ 
ture, spectroscopy, and scattering at low energies. It is be¬ 
ing learned quantitatively through lattice calculations of 
quark spin from the anomalous Ward identity [168|il00] 
that the smallness of the quark spin contribution in the 
nucleon is mainly due to the 17(1) anomaly, the same ori¬ 
gin which is responsible for the large rf mass. This cannot 
be understood with quark models without the chiral 17(1) 
anomaly. Similarly, relativistic quark models do not ex¬ 
plain th e larg e OAM obtained from the lattice calculation 
in Sec. 


5.1.2 Both the chiral quark model studies 


and lattice calculation of valence QCD ]1021I106| reveal 
that the level reversal of the positive and negative parity 
excited states of the nucleon, i.e. Pn(1440) (Roper res¬ 
onance) and S'n(1535), and the hyperfine splittings be¬ 
tween the decuplet and octet baryons are dominated by 
the meson-mediated flavor-spin interaction, not the gluon- 
mediated color-spin interaction. All of these point to the 
importance of the meson degree of freedom (qq pairs in the 
higher Fock space) which is missing in the quark model. 

To see how this comes about, one can follow Wilson’s 
renormalization group approach to effective theories Hi 1 1 
|141|. It is suggested bv Liu et al. una that the effective 
theory for baryons between the scale of 47r/ w , which is the 
scale of the meson size {Im ~ 0.2 fm), and ~ 300 MeV, 
which is the scale of a baryon size {Ib ~ 0.6 fm), should 
be a chiral quark model with renormalized couplings and 
renormalized meson, quark and gluon fields which pre¬ 
serve chiral symmetry. A schematic illustration for such 
division of scale:-[^]for QCD effective theories is illustrated 
in Fig. [3] 


QCD 

(V. V. 

.. 

0.2 fm (Ajt= 4ltf n) 

Chiral Effective Theory 
(V L . VL' A p. "-P- 3 !’ ■■■) 


0.6 fm 


Q b ~0.6 fm 


Chiral Perturbation Theory 

(? B . V B - n ' P' a c ■■■) 


Fig. 3. A schematic illustration of the the two-scale delineation 
of the effective theories. The shaded bars mark the positions of 
the cutoff scales Im and lg separating different effective theo¬ 
ries. 


This suggestion is based on the observation that mesons 
and baryon form factor assume a monopole and dipole 
form, respectively. Since the ttNN form factor is much 
softer than the pint form factor, it is suggested that the 
confinement scale of quarks in the baryon Ib is larger than 
Im ~ the confinement scale between the quark and an¬ 
tiquark in the meson. This is consistent with the large- 
N c approach to hadrons where the mesons are treated as 
point-like fields and the baryons emerge as solitons with 
a size of order unity in N c (165] . Taking Im from the 
p7T7r form factor gives Im ~ 0.2 fm. This is very close 
to the chiral symmetry breaking scale set by A x = 
Considering them to be the same, operators of low-lying 
meson fields become relevant operators below A x . As for 
the baryon confinement scale, Liu et al. take it to be the 
size characterizing the meson-baryon-baryon form factors. 
Defining the latters from the respective meson poles in the 
nucleon pseudoscalar, vector, and axial-vector form factors 
in a lattice calculation 11021 (see Fig. 17 in the reference), 
Liu et al. obtained Ib ~ 0.6 — 0.7 fm, satisfying Ib > Im- 
Thus, in between these two scales Im and Is , one could 
have coexistence of mesons and quarks in an effective the¬ 
ory for baryons. 

An outline is given by Liu et al. IM] to show how 
to construct a chiral effective theory for baryons. In the 
intermediate length scale between Im and Ib , one needs 
to separate the fermion and gauge fields into long-range 
ones and short-range ones 

ip =ipL+tps, A^ = Ag + Ag, ( 68 ) 

where V’l/V’S and A^/Ag represent the infrared/ultraviolet 
part of the quark and gauge fields, respectively, with mo¬ 
mentum components below/above 1/Im or A x . One adds 
irrelevant higher-dimensional operators to the ordinary 
QCD Lagrangian with coupling between bilinear quark 


11 We should point out that although two scales are adopted 
here, they are distinct from those of Manohar and Georgi ma¬ 
in the latter, the a - quark model does not make a distinction 
between the quark fields in the baryons and mesons. As such, 


there is an ambiguity of double counting of mesons and qq 
states. By making the quark-quark confinement length scale 
Ib larger than the quark-antiquark confinement length scale 
Im, one does not have this ambiguity. 
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fields and auxiliary fields as given in Ref. [98], interpret¬ 
ing these quark fields as the short-range ones, i.e. ips and 
ip s . Following the procedure by Li in Ref. [98] . one can in¬ 
tegrate out the short-range fields and perform the deriva¬ 
tive expansion to bosonize ips and ips- This leads to the 
Lagrangian with the following generic form: 


C-xQCD = C-QCD'iipLi^LiA^) + Cm^: Pi ai,G,- ■ ■) 

+ £<rq(.4’L, 'pL,K, O’, p,ai,G, ■ ■■). 

(69) 

Cqcd' includes the original form of QCD but in terms of 
the quark fields ip L ,ipL, and the long-range gauge field A £ 
with renormalized couplings. It also includes higher-order 
covariant derivatives [TSl] . Cm is the meson effective La¬ 
grangian, e.g. the one derived by Li jj33| which should in¬ 
clude the glueball field G. Finally, C aq gives the coupling 
between the ip L , ipL, G and mesons. As we see, in this 
intermediate scale, the quarks, gluons and mesons coexist 
and meson fields couple to the long-range quark fields. Go¬ 
ing further down below the baryon confinement scale 1 /Ib, 
one can integrate out ip L ■> V’i and A resulting in an ef¬ 
fective Lagrangian £(<F b,'Cb,'k , <j, p,a±,G, ■ ■ ■) in terms of 
the baryon and meson fields m■ This would correspond 
to an effective theory in the chiral perturbation theory. In 
order for the chiral symmetry to be preserved, the effective 
theory of baryons at the intermediate scale necessarily in¬ 
volves mesons in addition to the effective quark and gluon 
fields. This naturally leads to a chiral quark effective the¬ 
ory. 

Models like the little bag model with skyrmion outside 
the MIT bag [25], the cloudy bag model [151] and quark 
chiral soliton model [ 158] have the right degrees of freedom 
and, thus, could possibly delineate the pattern of division 
of the proton spin with large quark OAM contribution. In 
particular, the fact that the u and d OAM in the MIT bag 
and to some extent the LFCQM in Table 1 have different 
signs from those of the lattice calculation may well be due 
to the lack of meson contributions as demanded by chiral 
symmetry in the effective theory of baryons. 


7 Spin-orbit correlation 


We have seen that parton OAM constitutes an important 
part of the proton spin puzzle. Lorce and Pasquini m 
showed that another interesting and independent infor¬ 
mation about the nucleon spin structure is encoded in the 
correlation between the parton OAM and its own polariza¬ 
tion. This quantity escaped attention because it does not 
contribute to the nucleon spin budget, but naturally ap¬ 
pears once one considers all correlations allowed by parity. 
Parton spin-orbit correlation is in some sense the parity 
partner of parton OAM ESI- It is therefore hardly sur¬ 
prising that most of the results and discussions we had 
about parton OAM find a counterpart for parton spin- 
orbit correlation. 


For example, the quark kinetic and canonical spin- 
orbit operators are defined as mm 


r q 

°kin 


= J d 3 rip^"f 5 (r x iD)ip , 
Qan = J d 3 r- 0 t 75 (r x iD pure )ip. 


(70) 


Like the average quark kinetic OAM, the average quark 
kinetic spin-orbit correlation can be expressed in terms of 
moments of twist-2 and twist-3 GPDs m 

(C'kin) = \ J d xxH q (x,0,0) - | [F?{0) 

= - J dx x[G\{x, 0 , 0 ) + 2G\(x, 0 , 0 )]. (71) 

The parton spin-orbit correlation can also be naturally 
expressed as a phase-space integral mm 

(C q w G ) = 

f (72) 

J dxd 2 k±d 2 b± (b± x k±) z p++{x, k±, b±; W), 


where p+ G (x, k±, b±; W) is the helicity phase-space distri¬ 
bution describing the difference between the distributions 
of quarks/gluons with polarization parallel and antiparal¬ 
lel to the longitudinal direction. In terms of the GTMDs, 
this relation reads [117H85l[TT5] 

(Ctf) = J dxd 2 k ± ^G q 1 f(x,0,k ± ,0 ± ;W), ( 73 ) 

where once again the shape of the Wilson line W deter¬ 
mines whether the spin-orbit correlation is of the kinetic 
or canonical type. 

Like in the OAM case, the integrand of Eq. ( [72[ itself 
can consistently be interpreted as the phase-space density 
of quark/gluon spin-orbit correlation 

(C* G )(x,k ± ,b ± ) = (b± x fcjjz k±,b±;W). (74) 

From parity and time-reversal symmetries, it follows that 
the helicity phase-space distribution of longitudinally po¬ 
larized quarks/gluons inside a longitudinally polarized nu¬ 
cleon can be decomposed into four contributions 

pf + = pf G + (b ± ■ k ± ) p q ’ G + (b ± X k ± ) z p q ’ G 

+ (b 1 xfe 1 ) z (b 1 .fe 1 )pf, 

where pf G = pf G (x, fej_, (b± ■ k±) 2 , b^; W). The func¬ 
tions pi and p 2 are related to the Fourier transform of 
the GTMD G 14 , whereas p 3 and p 4 are related to the 
Fourier transform of the GTMD G n. Once again, the fac¬ 
tor (b± ■ k±) indicates that p 2 and pi are naive T-odd, i.e. 
they change sign under -f-t Wc, and arise therefore 
from the imaginary part of the corresponding GTMDs. 
Note also that the only contribution to the spin-orbit cor¬ 
relation surviving integration over 6 j_ or k± is p 3 . For a 
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straight Wilson line WY one obtains the density of JXY 
quark and gluon OAM [ySHIHT] provided that one integrates 
at least over k± [IIB- 

Using the extraction of polarized quark distributions 
from combined inclusive and semi-inclusive DIS data ffiBj 
in very good agreement with lattice calculations [21], Lorce 
found the estimates (C£ in ) ~ —0.8 and (C^ in ) ~ —0.55, 
meaning that quark spin and kinetic OAM are in aver¬ 
age antiparallel. On the contrary, the canonical version 
of the quark spin-orbit correlation appears to be positive 
in light-font quark models m, illustrating the impor¬ 
tance of the quark-gluon interaction. Here also it would 
be very interesting to compute these quantities using the 
large momentum approach of lattice gauge theory. 
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